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ILLUSTRATIONS 


Figure  1.  Ray  paths  for  the  Ps  phase  relative  to  the  P  phase  for  a  layer  over  a  half 
space.  Tp  and  Ts  are  the  travel  times  of  the  P  and  S  phases  with  the  same  ray 
parameter  through  the  layer  respectively.  Th  is  the  travel  time  differential  in  the  half 
space  for  the  two  rays  assuming  a  planar  wave  front. 

Figure  2.  Seismic  section  of  synthetic  receiver  functions  computed  for  a  range  of  ray 
parameter  from  0.040  to  0.080  by  raytracing  through  a  model  with  a  40  km  thick  layer 
(Vp=6.0  km/s,  Vs=3.5  km/s)  over  a  half  space  {Vp=8.0  km/s,  Vs=4.6  km/s). 

Figure  3.  Ray  paths  for  the  P2pis  phase  relative  to  the  P  phase  for  a  layer  over  a  half 
space.  Tp  and  Ts  are  the  travel  times  of  the  P  and  S  phases  with  same  ray  parameter 
through  the  layer  respectively.  Th  Is  the  travel  time  differential  in  the  half  space  for  the 
two  rays  assuming  a  planar  wave  front. 

Figure  4.  Velocity  spectrum  stacks  produced  from  the  synthetic  receiver  functions 
depicted  in  Figure  2.  Ps  stacks  are  shown  on  the  left  and  Pzpis  stacks  are  on  the 
right. 

Figure  5.  Stacks  of  the  synthetic  receiver  functions  depicted  in  Figure  2.  On  top  is  a 
straight  stack  with  no  time  correction  applied.  The  bottom  three  receiver  function 
stacks  are  computed  after  applying  the  appropriate  normal  moveout  correction  for  Ps, 
P2pis,  and  Pip2s  respectively  (from  top  to  bottom). 

Figure  6.  Moveout  curves  at  various  depths  from  25  km  to  700  km  for  the  Ps  phase 
computed  for  the  PREM  velocity  structure.  The  solid  lines  are  moveout  curves 


computed  for  a  single  layer  with  equivalent  average  slowness  as  PREM  to  the 
respective  depths.  The  dashed  line  is  the  moveout  computed  for  a  curved  ray  path  by 
ray  tracing  through  PREM.  The  depths  associated  with  each  curve  are  printed  on  the 
left. 

Figure  7.  Moveout  curves  at  various  depths  from  25  km  to  700  km  for  the  P2pis  phase 
computed  for  the  PREM  velocity  structure.  The  solid  lines  are  moveout  curves 
computed  for  a  single  layer  with  equivalent  average  slowness  as  PREM  to  the 
respective  depths.  The  dashed  line  is  the  moveout  computed  for  a  curved  ray  path  by 
ray  tracing  through  PREM.  The  depths  associated  with  each  curve  are  printed  on  the 
left. 

Figure  8.  A  representative  display  of  the  80  low  pass  filtered  (corner  frequency  of  0.3 
Hz)  PREM  synthetic  receiver  functions  used  to  compute  the  VSS  in  later  figures.  The 
moveout  curve  associated  with  Moho  and  upper  mantle  discontinuities  are  labeled. 
The  dashed  lines  are  the  true  moveout  curves  produced  by  raytracing  through  PREM. 
The  solid  lines  are  the  moveout  curves  computed  by  replacing  PREM  with  a  layer  over 
a  that  would  result  in  a  vertical  travel  time  equivalent  to  that  of  PREM. 

Figure  9.  Spherical  earth  velocity  spectrum  stacks  computed  using  the  PREM 
synthetics  depicted  in  Figure  8.  The  plot  to  the  left  was  computed  for  Ps  phases  while 
the  one  on  the  right  was  computed  using  the  P2pis  moveout  curve. 

Figure  10.  The  VSSc  produced  from  the  PREM  synthetics  depicted  in  Figure  8.  The 
VSSc  on  the  top  was  computed  using  PREM  P  and  S  velocities  as  the  reference 
model.  The  reference  model  used  to  produce  the  two  VSSc  on  the  bottom  used  P 
velocities  from  PREM  (left)  and  K8  (right)  and  a  value  of  1 .825  for  Rv. 
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Figure  1 1 .  Velocity  ratio  spectrum  stacks  (^ySS)  computed  using  PREM  (left)  and  K8 
(right)  reference  models  for  the  synthetic  receiver  functions  generated  for  the  PREM 
model. 

Figure  12.  A  representative  sampling  of  the  113  low  pass  filtered  individual  receiver 
functions  computed  for  OBN.  The  only  clear  arrivals  without  stacking  are  the  Ps 
conversion  from  the  Moho  at  5  seconds  and,  in  a  few  receiver  functions  the  P2pis  and 
Pip2s  from  the  Moho  at  20  and  23  seconds  respectively. 

Figure  13.  Synthetic  (solid  line)  and  observed  (dashed  line)  stacked  receiver 
functions  computed  for  the  IRIS/IDA  Seismographic  station  at  Obninsk.  The  synthetics 
pertain  to  the  crustal  structure  models  depicted  on  the  left  of  the  respective  receiver 
functions.  The  low  pass  filtered  response  (stack  of  the  receiver  functions  given  in 
Figure  12)  is  given  in  the  top  frame;  the  broad  band  results  are  in  the  bottom  frame. 

Figure  14.  Velocity  spectrum  stacks  produced  from  the  OBN  receiver  functions  shown 
on  Figure  8.  The  plot  on  the  left  was  computed  for  Ps  phases  while  the  one  on  the 
right  was  computed  using  the  P2pis  moveout  curve. 

Figure  15.  The  VSSc  produced  from  113  OBN  receiver  functions  using  PREM  as  the 
reference  model. 

Figure  16.  Velocity  ratio  spectrum  stack  (R^^S)  from  the  OBN  receiver  functions 
computed  using  PREM  as  the  reference  model. 


Figure  17.  VSSc  computed  from  OBN  receiver  functions  similar  to  that  shown  in 
Figure  15,  except  in  this  case,  the  VSSc  was  computed  using  a  value  of  1.78  for  Rv. 

Figure  1 8.  Velocity  spectrum  stacks  produced  from  receiver  functions  computed  from 
seismograms  recorded  at  ARU.  The  plot  on  the  left  was  computed  for  Ps  phases 
while  the  one  on  the  right  was  computed  using  the  P2pis  moveout  curve. 

Figure  19.  Velocity  ratio  spectrum  stack  (RySS)  from  the  ARU  receiver  functions 
computed  using  PREM  as  the  reference  model. 

Figure  20.  VSSc  computed  from  ARU  receiver  functions  similar  to  that  shown  in 
Figure  19,  except  in  this  case,  the  VSSc  was  computed  using  a  value  of  1.78  for  Rv. 


SUMMARY 


In  order  to  improve  the  signal  to  noise  ratio  of  receiver  function  data,  it  is  typical  to 
stack  receiver  functions  calculated  from  events  at  similar  distances  and  back  azimuths. 
We  have  adapted  the  velocity  spectrum  stacking  (VSS)  technique,  used  extensively  in 
reflection  seismology,  to  the  receiver  function  method  in  order  to  stack  data  with 
different  ray  parameters,  thereby  improving  further  the  signal  to  noise  ratio.  Perhaps 
more  importantly,  by  producing  the  velocity  spectrum  stacks  we  take  advantage  of  the 
differences  in  the  shapes  of  the  moveout  curves  of  converted  phases  and 
reverberations  to  identify  and  separate  the  various  phases  and  to  infer  velocity 
structure.  Through  conventional  receiver  function  techniques  we  have  modeled  the 
crustal  structure  beneath  the  Russian  IRiS/lOA  seismographic  station  at  Obninsk, 
Russia  (OBN).  This  model  includes  a  2  km  thick  low  velocity  surface  layer  and  47  km 
depth  to  Moho  with  relatively  uncomplicated  crustal  structure.  By  computing  velocity 
spectrum  stacks  from  the  observed  seismograms  at  OBN  and  Arti  (ARU),  Russia,  and 
comparing  them  with  those  produced  from  PREM  synthetics  we  have  identified  Ps 
phases  from  the  400  and  670  km  discontinuities.  We  find  that  these  phases  can  be 
interpreted  satisfactorily  in  terms  of  a  mantle  structure  characterized  by  higher  upper 
mantle  S  velocities,  and  possibly  deeper  discontinuities  than  in  PREM  and  in  these 
data  we  find  no  evidence  of  a  200  km  discontinuity  beneath  either  of  these  stations. 
By  comparison  of  VSS  between  OBN  and  ARU  we  suspect  that  the  670  km 
discontinuity  is  about  2  km  shallower  beneath  ARU  than  OBN. 
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1.  INTRODUCTION 


A  commonly  used  technique  to  estimate  crust  and  upper  mantle  structure  from  a 
single  three-component  seismographic  station  is  to  compute  and  interpret  “receiver 
functions”  (e.g.,  Langston  1989,  1981,  1979,  1977;  Owens,  et  al.  1984,  1987,  1988; 
Owens  and  Crosson  1988),  wherein  the  horizontal  components  are  deconvolved  by 
the  vertical  component  to  produce  a  trace  consisting  primarily  of  Ps  conversions  and 
converted  S-wave  reverberations.  The  technique  has  been  successfully  extended  to 
arrays  of  broadband  portable  stations  by  Owens  et  al.  (1988a,b).  To  improve  the 
signal  to  noise  ratio,  receiver  functions  can  be  binned  by  ray  parameter  and  back 
azimuth  and  stacked  (Owens  et  al.  1983;  Owens  1984).  In  areas  with  flat  geological 
structure  the  receiver  functions  show  little  or  no  azimuthal  dependence  and  can  be 
stacked  at  common  ray  parameter  for  all  azimuths.  However,  if  we  wish  to  stack  traces 
with  different  ray  parameters,  we  must  first  correct  for  the  differences  in  relative  arrival 
times. 

The  “velocity  spectrum  stack”  (VSS)  is  a  useful  tool  for  stacking  reflection  data 
within  a  range  of  ray  parameters  in  multichannel  studies  (e.g.  Yilmaz  1987).  The 
functional  dependence  of  arrival  times  on  ray  parameter  p,  relative  to  a  reference 
phase  with  ray  parameter  po  is  called  the  “moveout”.  The  “normal  moveout  correction” 
(NMO)  then  refers  to  the  time  adjustment  necessary  to  correct  the  arrival  time  to  what 
would  have  been  observed  from  a  vertically  incident  ray,  irrespective  of  amplitude, 
assuming  a  given  velocity  structure.  The  “velocity  spectrum  stack”  (VSS)  is  a  contour 
map  of  amplitudes  across  constant  velocity  stacks  (produced  by  stacking  the  observed 
records  after  NMO  assuming  a  uniform  velocity)  in  the  velocity-time  plane  (e.g.  Sheriff 
1982).  A  phase  present  in  the  receiver  functions  is  thus  enhanced  in  the  VSS  if  the 
appropriate  NMO  correction  is  made.  The  enhancement  will  be  most  effective  for  a 
value  of  velocity  matching  the  “true”  mean  velocity  sampled  by  the  phase,  it  is  most 
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appropriate  at  this  point  to  think  of  the  velocity  structure  as  a  function  of  time  since 
arrival  time  is  observed  whereas  depth  will  be  computed  after  the  velocity  structure  is 
determined  (Yilmaz  19b/).  Because  of  differences  in  the  shapes  of  their  moveout 
curves,  separate  stacks  must  be  produced  for  each  of  the  prominent  phases  present  in 
the  receiver  functions.  Therefore,  the  velocity  spectrum  stacks  can  be  used  to 
disc  «guish  between  phases  as  well  as  to  infer  velocity  structure. 

The  examples  in  the  following  sections,  describing  the  production  of  VSS,  use 
radial  components  of  the  receiver  functions.  Since  the  production  of  the  VSS  depends 
only  on  travel  time,  the  tangential  component  VSS  can  be  produced  using  the  same 
procedure,  however  for  truly  one-dimensional  structure  there  will  be  no  energy  on  the 
tangential  components.  Since  the  emphasis  of  this  paper  is  the  technique  most  of  the 
explanations  will  use  synthetic  data.  We  will,  however,  discuss  VSS  produced  from 
three  component  broad  band  seismograms  recorded  at  Russia  IRiS/IDA  stations  at 
Obninsk  (OBN)  and  Arti  (ARU). 

2.  METHOD 

Figure  1  illustrates  the  geometry  of  the  most  significant  type  of  phase  observed  in 
receiver  function  studies  —  the  P  to  S  conversion  (Ps)  generated  when  the  wave 
crosses  an  interface  —  for  a  layer  over  a  half  space.  We  assume  a  planar  incoming 
wave  front  (a  typical  assumption  in  receiver  function  studies)  in  deriving  the  following 
equations.  The  time  delay  for  the  Ps  arrival  relative  to  that  of  the  P  arrival  ATps(p)  is 
given  by: 

ATp.(p)  =  Ts+Th*Tp  (1) 

ATp,(2,p,Vs,Vp)  =  z (Vv^-p2  -  VV-,?.p2)  (2) 
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In  the  above  equations:  Ts,  Th  and  Tp  are  travel  times  along  the  paths  labeled  in 
Figure  1 ;  Vs  and  Vp  are  the  average  S  and  P  velocities  in  the  layer,  respectively;  p  is 
the  ray  parameter;  and  z  is  the  depth  to  the  interface.  In  terms  of  the  vertical  travel  time 
of  Ps  relative  to  P  (ATpso)  through  the  layer  and  the  velocity  ratio  Rv  =  Vp  /  Vg,  we 
have: 


ATps(ATpso.P.Vs.flv^  =  %AlPa0(Vl-p2V|  -  )  (3) 

Hv  -  1 

Note  that  this  equation  depends  only  on  ATpso.  p.  Vs  and  an  assumed  value  for  Rv 
(e.g.  Rv  =V3  for  a  Poisson  solid). 

Figure  2  depicts  a  set  of  synthetic  receiver  functions  generated  for  a  layer  (Vp  =  6.0 
km/s,  Vs  =  3.5  km/s)  over  a  half  space  (Vp  =  8.0  km/s,  Vs  =  4.6  km/s,  and  z  *  40  km)  for 
a  range  of  ray  parameters.  We  see  that  the  Ps  phase  is  delayed  with  increasing  ray 
parameter  relative  to  the  initial  P-wave,  illustrating  the  Ps  moveout.  All  phases 
(reverberations)  following  Ps  are  advanced  in  arrival  time  relative  to  the  P  phase  with 
increasing  ray  parameter.  The  next  prominent  phase  after  the  Ps  arrival  is  composed 
of  the  sum  of  all  reverberations  in  which  there  are  two  P  legs  and  one  S  leg.  Figure  3 
depicts  the  travel  path  of  the  Ppps  phase  in  which  the  first  two  branches  are  P  waves 
and  the  final  branch  is  an  S.  For  near  vertical  incidence  the  reverberations  ending  in 
an  S  leg  will  contribute  much  more  energy  to  the  horizontal  components  than  those 
ending  in  a  P  (i.e.  Pspp,  and  Ppsp).  For  convenience  of  notation  we  will  adopt  a 
naming  convention  for  the  sum  of  these  reverberations  as  Pnpms,  in  which  n  is  the 
number  of  p  legs  and  m  is  the  number  of  s  legs  (if  n  or  m  is  zero  it  is  omitted).  The 
phase  described  will  be  referred  to  as  Pzpis.  The  next  phase  on  Figure  2,  Pip2s, 
composed  of  the  sum  of  reverberations  with  two  S  legs  and  one  P,  has  two 
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contributions  with  final  S  legs  (the  Psps  and  Ppss)  and  a  less  significant  Pssp 
contribution. 

We  take  advantage  of  this  difference  in  the  shape  of  the  moveout  curve  to 
distinguish  reverberations  from  the  Ps  phase.  Equation  4  gives  the  time  delay 
(ATp2pis)  for  the  P2pis  phase  relative  to  the  P  arrival  for  a  layer  over  a  half  space: 


ATp2p,s(ATp2p,so.p.Vs.flv)  =  +  Vflv-^-p2V| )  (4) 

In  like  fashion,  we  can  derive  moveout  equations  for  the  Pap,  Pip2s  and  Pas 
phases.  However  Pap  and  Pas  have  small  amplitudes  (even  after  stacking  a  very 
large  number  of  events)  so  are  usually  of  tittle  significance  in  interpretations.  Pip2s, 
on  the  other  hand,  has  reversed  polarity  and  is  easily  distinguishable  from  Ps  and 
P2pis,  so  we  found  it  a  useful  phase  in  the  interpretation  of  receiver  structure  at  OBN. 
The  moveout  for  this  phase  is  given  by: 


ATp,pa(ATp,pao.P.Vs.ny)  -  ^  Vi  -  p»V|  (5) 

Constant  velocity  stacks  are  produced  by  averaging  along  the  moveout  curve  the 
amplitudes  of  N  receiver  functions  with  various  ray  parameters. 

S(AT4K).V,)  =  1  X  fi(AT«(AToo.Pi.Vs.«v^)  (6) 

i.i 

In  this  equation,  <I>  is  the  type  of  phase  (e.g.  Ps  or  P2pis);  S(AT<i>o  Vg)  is  the 
averaged  amplitude  at  a  given  zero  offset  time  and  S*wave  velocity; 
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fi{AT<i,{AT«i,o.Pi.Vs.Wi')}  is  the  amplitude  of  the  |th  trace  at  the  computed  moveout  time. 
AT<t(AT4io,Pi,Vs,/?v'),  for  a  given  wave  type  (^).  If  the  moveout  time  falls  between  two 
samples  we  linearly  interpolate  a  value  for  fj{AT<[i(AT4>o>Pi.Vs,/?v)}-  After  producing 
constant  velocity  stacks  for  the  range  of  all  reasonable  velocities,  we  contour  the 
amplitudes  in  the  velocity*time  plane  to  produce  the  VSS.  The  Ps  conversion  or  P2pis 
reverberation  on  their  respective  velocity  spectrum  stacks  will  appear  as  positive 
ridges  (negative  for  a  velocity  inversion)  elongated  parallel  to  the  velocity  axis.  The 
velocity  structure  beneath  a  station  can  then  be  inferred  by  selecting  the  time  and 
velocity  of  the  highest  amplitude  on  each  ridge. 

An  alternative  to  the  summation  of  receiver  function  amplitudes  along  the  moveout 
curve  to  generate  a  VSS  is  the  normalized  summation  of  the  zero  lag  crosscorrelation 
of  all  the  receiver  functions  in  a  small  window  (usually  the  wavelength  of  the  expected 
arrival)  centered  about  the  moveout  curve  (Yilmaz,  1987).  We  have  found  that  the 
advantages  gained  by  this  approach  are  usually  cosmetic  and  occasional  instabilities 
arise  when  it  Is  applied  to  noise  free  synthetics  or  data  in  which  the  moveout  curves  of 
the  Ps  phase  and  reverberations  cross.  Although  in  some  cases  it  may  be  beneficial 
to  produce  VSS  by  crosscorrelation,  we  feel  that  this  does  not  add  substantially  to  the 
major  points  of  this  paper  and  refer  the  reader  to  discussions  in  reflection  seismology 
texts  (e.g.  Yilmaz,  1987). 

3.  SYNTHETIC  EXAMPLE 

Figure  4  depicts  VSS  produced  for  the  Ps  and  Pzpis  phases  using  the  synthetic 
receiver  functions  shown  in  Figure  2.  Upon  inspection  of  the  Ps  stack  (left)  we  observe 
good  time  resolution  for  Ps  near  5  seconds  but  poor  velocity  resolution.  This  phase 
also  appears  on  the  Pzpis  stack  (right),  but  the  peak  is  not  as  sharp  and  does  not 
have  as  large  an  amplitude  as  on  the  Ps  stack.  Pzpis  (at  18  s)  is  only  observed  on  the 
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P2pis  Stack  and  has  much  better  velocity  resolution  than  Ps.  It  is  not  surprising  that 
we  are  able  to  pick  Vs=3.5  km/s  (the  velocity  used  to  compute  the  synthetics  in  Figure 
2)  more  accurately  from  the  P2pis  stack,  sin(»  Figure  2  shows  twice  as  much  moveout 
for  P2P  IS  than  for  Ps.  We  can  use  this  velocity  together  with  the  approximate  4.8  s 
arrival  time  on  the  Ps  stack  to  compute  the  thickness  of  the  layer.  The  arrival  times  on 
the  VSS  are  for  vertically  incident  rays  (p=0),  so  equation  (2)  becomes: 

z  =  Tps/{1/Vs-1/Vp}  (7) 

For  this  shallow  structure,  we  assume  a  poisson  solid  {R\/  =V3)  and  solving  for  z  we 
confirm  the  40  km  depth  to  the  interface. 

Figure  5  depicts  single  stacks  of  the  receiver  functions  shown  in  Figure  2.  The  top 
stack  has  no  moveout  applied  -  the  next  three  are  stacked  using  the  respective  Ps, 
P2pis,  and  Pip2s  moveout  curves  (from  top  to  bottom).  In  each  case  the  respective 
velocity  depth  (time)  function  was  picked  by  observation  of  the  stacking  amplitudes  on 
the  corresponding  VSS  (Figure  4).  The  P,  Ps.  P2pis,  and  Pip2s  arrivals  are  at  5,  10. 
23  and  28  seconds  respectively.  Figure  5  clearly  illustrates  the  fact  that  the  various 
arrivals  are  substantially  enhanced  when  stacked  along  the  appropriate  moveout 
curve.  An  added  bonus  is  the  annihilation  of  the  reverberations  (Pzpis  and  Pip2s)  on 
the  Ps  stack;  conversely  the  Ps  phase  is  greatly  diminished  on  the  two  reverberation 
stacks.  We  conclude  that  by  producing  the  stacks  with  normal  moveout  we  may 
identify  arrivals  that  would  otherwise  not  be  observed' and,  in  the  process,  avoid 
mislabeling  other  phases. 

Figures  6  arid  7  depict  the  shapes  of  the  Ps  and  P2pis  moveout  curves  computed 
for  various  depths  using  the  PREM  velocity  model  (Dziewonski  and  Anderson,  1981) 
over  the  range  of  ray  parameters  typical  of  teleseismic  P-arrivals  used  in  receiver 
function  studies  (0.04  to  0.08  s/km  for  epicentral  distances  of  30°  to  90°).  For  each 
given  depth,  the  solid  curve  was  computed  for  a  single  layer  with  average  slowness  (P 
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and  S)  equivalent  to  that  of  PREM.  The  dashed  lines  are  the  moveout  curves 
computed  to  include  the  additional  time  delay  resulting  from  the  curvature  of  the  ray 
path  through  PREM,  discussion  of  which  we  defer  to  a  later  section.  At  this  time,  we 
shall  merely  point  out  that  the  curvature  of  the  ray  paths  affects  significantly  the  shape 
of  the  moveout  curves  only  for  interface  depths  greater  than  about  200  km. 

For  an  interface  depth  of  50  km,  there  are  0.75  s  of  moveout  over  the  teleseismic 
range  of  ray  parameter.  There  are  about  2.5  s  of  moveout  calculated  over  the  same 
range  of  ray  parameter  for  the  P2pis  phase  at  50  km  depth.  In  Figure  4  we  observed 
that  for  an  even  shallower  depth  (40  km),  with  less  moveout,  reasonable  signals 
appeared  on  the  VSS  for  both  phases.  For  the  Ps  phase,  there  are  1 .8  s  of  moveout  at 
100  km  over  the  range  of  ray  parameters,  and  as  much  as  16.1  s  at  600  km.  It  is  clear 
that  with  a  finite  number  of  traces  to  stack,  VSS  produced  for  the  Ps  phase  will  be 
much  more  useful  in  the  interpretation  of  upper  mantle  structure  than  of  crustal 
structure.  On  the  other  hand,  about  1 .25  seconds  of  moveout  is  expected  for  the  P2pis 
over  the  given  range  of  ray  parameter  at  25  km  depth.  Because  of  the  greater  amount 
of  moveout  calculated  for  this  phase,  its  VSS  may  prove  to  be  more  valuable  in  the 
interpretation  of  shallow  structure.  On  the  other  hand,  the  longer  ray  paths  of  the 
reverberations  (as  opposed  tQ  the  direct  Ps)  make  such  phases  more  sensitive  to 
heterogeneities  in  the  near  surface  structure.  As  a  result  interpretation  of  these 
phases  may  often  be  more  difficult  (e.g.,  Owens  et  al.  1984). 

It  is  clear  from  Figure  6  that  to  produce  reliable  VSS  for  Ps  phases  from  shallow 
layers  (100  km  or  less),  data  from  the  full  range  of  ray  parameters  are  necessary.  For 
the  Ps  phases  from  deeper  interfaces  and  the  P2pis  phase  at  ail  depths,  it  appears 
that  there  is  sufficient  moveout  along  the  curves  to  produce  VSS  from  data  with  poorer 
distribution  in  slowness.  In  cases  where  there  is  not  enough  distribution  of  data  to  use 
the  VSS  method  to  estimate  velocity  structure,  coherence  between  individual  receiver 
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functions  may  still  be  improved,  prior  to  stackir^,  by  making  a  moveout  correction  (not 
necessarily  normal  moveout)  using  velocities  from  a  regional  model. 

4.  EARTH  FLATTENING  TRANSFORMATION 

Traditionally,  receiver  functions  are  used  to  model  crustal  structure  and  it  is 
adequate  to  assume  a  flat  earth.  In  the  case  of  upper  mantle  discontinuities,  the 
curvature  of  the  earth  becomes  a  significant  factor.  Because  all  the  moveout 
equations  derived  above  were  based  on  a  fiat-layered  earth,  we  may  use  a  classical 
earth  flattening  transformation  (EFT.  Biswas  and  Knopoff  1970;  Aki  and  Richards, 
1980)  to  map  a  flat  earth  structure  into  the  spherical  earth  structure  that  will  yield  the 
equivalent  seismograms.  The  transformation  is  not  exact  for  P-SV  waveforms,  but  it  is 
exact  for  their  travel  times  (Chapman  1973).  Because  the  purpose  of  VSS  is  to  sum 
the  amplitudes  of  the  receiver  functions  along  the  moveout  curves  in  order  to  find  the 
velocity  and  time  delay  for  which  the  pulse  of  any  given  shape  stacks  most  coherently, 
travel  time  preserving  transformations  will  not  reduce  significantly  the  effectiveness  of 
the  VSS  method. 

The  earth  flattening  transformation  is  applied  to  each  grid  cell  of  the  VSS  in  order 
to  convert  the  flat  earth  velocities  (Vf)  to  spherical  earth  velocities  (Vsph)  according  to; 

Vf  =  Vsph(r/re)  .  (8) 

where  re  is  the  radius  of  the  earth.  Here,  r  is  the  radius  corresponding  to  the  grid  cell 
and  is  related  to  the  flat  earth  depth  (z)  by:  > 


r  =  r,  e^« , 


(9) 
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where  z,  for  each  grid  cell  of  the  time  versus  velocity  VSS,  is  determined  by  equation 
(7)  for  Ps  velocity  spectrum  stacks.  For  the  P2pis  phase  z  is  determined  by: 

z  =  Tps/(1/Vs+1/Vp)  (10) 

By  performing  this  transformation  on  a  flat  earth  VSS  (VSSf)  generated  similarly  to  that 
of  Figure  4  we  obtain  the  spherical  earth  VSS  (VSSsph)- 

Figure  8  shows  synthetic  receiver  responses  of  a  planar  P  wave  front  passing 
through  PREM  over  the  range  of  ray  parameter  from  0.04  to  0.08  s/km.  Each  of  these 
receiver  functions  was  low  pass  filtered  with  a  corner  frequency  of  0.3  Hz  in  order  to 
match  the  observations  at  OBN  discussed  below.  Figure  9  depicts  the  VSSsph 
computed  for  these  synthetics.  We  used  a  uniform  Rv  of  1.825  in  computing  the 
VSSsph  which  is  the  appropriate  average  value  for  PREM  to  a  depth  of  670  km.  The 
effect  of  this  approximation — as  opposed  to  using  a  depth-dependent  ratio —  will  be 
the  same  on  VSSsph  computed  for  synthetics  and  for  those  computed  for  observed 
data.  Letting  this  ratio  vary  with  depth  would  add  another  free  parameter  and 
unnecessarily  complicate  the  discussion.  The  shallower  arrivals  (for  which  a  lower 
value  of  R\^  is  appropriate)  will  stack  at  a  slightly  lower-than-expected  velocity  as  a 
result  of  this  approximation.  The  contour  intervals  for  Figure  9  were  chosen  to 
emphasize  the  upper  mantle  discontinuities,  so  that  crustal  phases  arriving  in  the  first 
25  seconds  of  the  plot  are  saturated.  The  Ps  conversions  from  the  400  and  670  km 
discontinuities  appear  at  about  41  seconds  and  63  seconds  respectively  on  the  Ps 
velocity  spectrum  (left).  Because  the  Ps  phase  from  the  200  km  discontinuity  is  on  the 
trailing  edge  of  the  Moho  reverberations  at  20  seconds  there  is  no  distinct  arrival  for  it 
on  the  VSS;  however  the  corresponding  Pzpis  phase  appears  strongly  at  76  seconds 
on  the  P2pis  VSSsph- 
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5.  EFFECTS  OF  RAY  PATHS  ON  VSS  WHEN  APPLIED  TO  UPPER 
MANTLE  DISCONTINUITIES 


In  the  previous  section,  we  outlined  the  correction  necessary  to  modify  the  VSS 
technique  for  application  to  deep  structure.  In  this  section  we  outline  the  procedures  to 
correct  for  errors  introduced  resulting  from  the  following  'ray  path'  assumptions;  1) 
typically  in  VSS  studies  the  effect  of  curvature  of  ray  path  is  ignored;  2)  in  both  receiver 
function  and  VSS  studies  a  planar  incoming  wave  front  is  assumed. 

In  the  process  of  mapping  the  VSSf  into  the  VSSsph  we  computed  the  radius  for 
every  grid  cell  on  the  VSS,  so  that  we  can  plot  the  result  as  a  velocity  spectrum  stack 
in  the  spherical-earth-depth  versus  velocity  domain  (VSSz)-  However,  returning  to 
Figures  6  and  7,  we  observe  that  the  true  moveout  curves”  produced  by  ray  tracing 
through  PREM  exhibit  a  greater  amount  of  moveout  than  those  computed  by  replacing 
the  PREM  structure  above  the  interface  with  a  single  layer  of  equal  average  vertical 
slowness.  Because  the  amount  of  moveout  is  directly  proportional  to  the  increase  in 
velocity,  it  is  clear  that  computing  moveout  using  the  single  layer  assumption  outlined 
above  will  result  in  a  high  estimate  of  velocity  when  applied  to  data  from  a  layered 
structure.  We  notice  in  Figure  6  and  7  that  the  true  moveout  curves"  and  the  single 
layer  moveout  curves  merge  to  the  same  zero  offset  time  delay.  As  a  result  the  time 
delays  observed  on  the  VSSsph  (Figure  9)  are  true  zero  offset  travel  times  but  the 
velocities  are  overestimated,  which  would  result  in  an  overestimate  of  the  depths  to  the 
upper  mantle  discontinuities.  This  should  affect  observed  data  and  synthetics  in  the 
same  way,  so  that  the  procedure  should  be  adequate  to  compare  synthetics  and 
observed  data;  it  will  of  course  work  well  for  shallow  layers.  To  estimate  more 
accurately  the  velocity  structure  directly  from  VSS,  we  have  refined  our  method  of 
producing  VSSz  by  applying  normal  moveout  corrections  computed  for  curved  ray 
paths  (VSSc). 


The  curved  ray  path  normal  moveout  correction  can  be  computed  by  ray  tracing 
through  a  reference  model.  Because  we  do  not  need  the  amplitude  information 
provided  by  ray  tracing,  it  is  quicker  to  compute  moveout  corrections  by  integrating  the 
relative  time  delay  (equations  2  through  4)  due  to  each  infinitesimally  thin  layer  (dz)  to 
the  given  depth  (zi). 


ATp8(z,p,Vs,Vp)  = 


2)-p2-Vvi?(z)'^-p2  )dz 


(11) 


Equation  (11)  is  used  to  compute  the  ATo  values  required  by  equation  (6)  to  produce 
“reference  velocity  model  stacks"  (RVMS)  in  the  same  manner  that  constant  velocity 
stacks  were  produced  in  previous  sections.  By  multiplying  all  velocities  of  the 
reference  model  by  a  constant  fraction  we  compute  a  "fractional  reference  model".  A 
VSSc  is  then  constructed  by  contouring  RVMS  produced  for  a  range  of  different 
fractional  reference  models.  The  axes  of  the  VSSc  are  depth  versus  "fraction  of 
reference  model",  therefore  a  phase  identified  on  the  VSSc  has  a  velocity  which  we 
describe  as  a  certain  fraction  or  percentage  above  or  below  the  reference  model. 

The  assumption  of  a  planar  wave  front  (i.e.  that  the  P  and  the  Ps  phases  have  the 
same  ray  parameter)  results  in  a  slight  error  in  the  estimate  of  the  shape  of  the 
moveout  curve.  This  error  behaves  similarly  to  that  described  above  for  the  curved  ray 
path  correction  but  results  in  only  about  a  fifth  as  much  error  in  the  estimated  velocity, 
and  therefore  would  cause  the  estimate  for  the  depth  to  the  670  km  discontinuity  to  be 
shallow  by  about  5  km.  For  shallower  layers,  the  error  in  depth  estimates  resulting 
from  this  assumption  is  negligible.  To  image  accurately  the  670  km  discontinuity  we 
must  compute  the  moveout  curves  by  ray  tracing  through  reference  models.  In 
constructing  the  VSSj  by  this  technique  we  found  that  for  certain  fractional  reference 
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models,  the  rays  with  large  ray  parameters  would  turn  before  reaching  the  surface. 
This  information  alone  can  be  used  to  exclude  certain  fractional  reference  models  from 
the  set  of  possible  solutions. 

Figure  10  shows  a  Ps  VSSc  for  the  same  PREM  synthetics  used  in  computing 
Figure  9.  The  top  VSSc  was  computed  using  the  PREM  P  and  S  velocities.  The  other 
two  plots  (bottom  left  and  right,  respectively)  were  computed  using  the  PREM  and  the 
K8  P'Velocity  model  for  northwestern  Eurasia  (Given  and  Helmberger,  1980)  and  a 
value  of  Rv-  1.825.  The  fact  that  the  leftmost  and  middle  plots  are  virtually  identical 
encourages  us  to  believe  that  a  model  such  as  K8  would  be  a  reasonable  reference 
model  provided  that  a  reasonable  value  for  is  used.  The  peaks  associated  with 
each  arrival  on  the  two  PREM  reference  model  VSSc  (left  and  center)  line  up  along  the 
value  of  1.0  on  the  horizontal  axis  and  appear  at  the  expected  670  km  depth. 
However  this  phase  appears  at  about  a  fractional  velocity  model  value  of  0.99  and  at  a 
slightly  greater  depth  on  the  K8  VSSc.  The  average  PREM  velocity  computed  to  this 
depth  is  about  .985  that  of  K8,  which  is  within  the  peak  observed  in  K8  reference 
model  VSSc  of  Figure  10.  This  small  difference  in  velocity  results  in  the  slight 
discrepancy  in  depth  when  applied  to  the  large  time  delays  associated  with  upper 
mantle  discontinuities. 

By  following  the  steps  outlined  above  to  generate  the  VSSc  except  holding  the 
model  constant  and  varying  the  velocity  ratio  (Ry ),  we  can  produce  a  spectrum  of 
stacked  receiver  functions  (R^SS)  which  will  be  most  coherent  at  the  proper  value  of 
Rv .  Figure  1 1  depicts  the  R^SS  computed  from  the  above  PREM  synthetics  using 
PREM  (left)  and  K8  (right)  as  reference  models.  For  this  synthetic  example,  we  find 
indeed  that  for  both  reference  models  the  expected  value  of  R^  =  1.82  to  1.825 
produces  the  greatest  amplitude  on  the  RySS. 


12 


6.  CRUSTAL  STRUCTURE  AT  OBN 


We  computed  Velocity  Spectrum  Stacks  for  data  recorded  at  the  Russian  station  at 
Obninsk  (OBN).  We  use  conventional  receiver  response  interpretation  to  constrain  the 
crustal  structure  and  use  the  VSS  method  in  the  following  section  to  look  at  upper 
mantle  discontinuities. 

Receiver  functions  for  the  IRIS/IDA  station  at  OBN  were  computed  using  data 
collected  in  1989-90  (Gurrola,  et  al.,  1990a,b).  The  station  is  equipped  with  a 
broadband  three-component  system  with  response  nominally  flat  with  respect  to 
velocity  from  approximately  3  mHz  to  5  Hz.  We  used  teleseismic  P  and  PP  phases 
which,  due  to  the  uneven  distribution  of  source  regions  during  the  one  year  period 
covered  by  the  data,  primarily  sample  the  northeast  and  southeast  quadrants.  We 
found  it  useful  to  high-pass  filter  the  seismograms  in  order  to  counter  the  effects  of 
occasional  nonlinear  noise  problems  at  frequencies  lower  than  20  mHz. 

The  broadband  OBN  receiver  functions  are  dominated  by  reverberations  within  a 
shallow  surface  layer.  In  order  to  identify  phases  from  deeper  layers,  we  have 
reduced  the  contribution  of  the  near  surface  layer  by  low-pass  filtering  these  data  with 
a  phaseless  Gaussian  filter  (with  a  half  power  width  of  0.6  Hz,  Figure  12).  The  velocity 
structure  was  determined  using  the  smooth  inversion  of  Ammon  (Ammon,  et  al.,  1990), 
which  employs  the  reflection  matrix  method  (Kennett,  1983;  and  Randall,  1989).  The 
simplest  model  that  we  could  construct  which  satisfies  both  the  broadband  and  the 
high  frequency  data  includes  a  low  velocity  surface  layer  of  no  more  than  2.5  km 
thickness  and  a  rather  smoothly  increasing  velocity  gradient  to  the  47  km  deep  Moho 
as  shown  on  Figure  13  (Gurrola  et  al.  1990a). 


13 


7.  UPPER  MANTLE  DISCONTINUITIES  AT  OBN 


VSSsph  produced  from  OBN  receiver  functions  (Figure  14)  exhibit  clear  arrivals 
from  the  upper  mantle  discontinuities.  The  Ps  VSSsph  is  depicted  on  the  left  and  the 
P2pis  VSSsph  is  on  the  right  of  the  figure.  We  observe  Ps  and  P2pis  phases  from  the 
Moho  at  about  5  and  20  seconds  respectively.  The  contour  interval  was  chosen  to 
illustrate  best  the  upper  mantle  arrivals  not  observable  on  the  individual  receiver 
functions,  so  that  Moho  arrivals  are  not  well  defined  on  these  plots.  We  used  the  same 
value  for  Rv  (1 .825)  as  used  in  the  PREM  VSS  discussed  above. 

The  observed  Ps  phase  from  the  400  km  discontinuity  (at  41  seconds  on  the 
VSSsph  of  Figure  14)  is  larger  in  amplitude  than  the  Ps  phase  from  the  670  km 
discontinuity.  This  is  the  opposite  of  our  observations  of  the  PREM  synthetics  in  which 
the  Ps  from  the  400  km  discontinuity  was  2/3  the  amplitude  of  that  of  the  670  km 
discontinuity  (Figure  9).  This  result  is  consistent  with  the  larger  velocity  contrast  for  the 
400  km  discontinuity  suggested  by  the  K8  continental  model  of  Given  and  Helmberger 
(1980).  Because  K8  is  only  specified  in  terms  of  P  velocities,  synthetics  produced  for  it 
will  be  heavily  dependent  upon  the  necessary  assumption  of  S  velocity  and  density, 
therefore  we  did  not  produce  a  synthetic  VSSsph  ^or  this  model.  The  Ps  arrival  from 
the  670  km  discontinuity  is  similar  in  amplitude  on  both  the  OBN  stacks  and  the  PREM 
synthetics,  which  implies  a  similar  velocity  contrast.  This  phase  appears  to  arrive 
slightly  earlier  in  time  and  at  a  higher  velocity  in  the  OBN  VSSsph  than  on  the  PREM 
VSSsph-  The  time  delay  between  the  400  and  670  km  dfscontinuities  is  smaller  than 
observed  in  PREM,  which  is  consistent  with  a  conclusion  reached  by  Vinnik,  based  on 
independent  analysis  of  OBN  data  (L.  Vinnik,  1990,  personal  comm.).  From  Figure  15 
(the  VSSc  computed  for  the  Ps  moveout  curve)  we  observe  that  the  expected  Ps 
phase  from  the  670  km  discontinuity  actually  stacks  best  as  a  phase  from  a  depth  of 
665  km  with  a  velocity  of  about  1 .04  times  that  calculated  from  PREM.  From  this  we 
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may  infer  that  this  discontinuity  is  either  5  km  shallower  beneath  OBN  than  predicted 
by  PREM  or  that  the  value  of  Rv  beneath  OBN  is  lower  than  in  PREM.  Figure  16 
depicts  RvSS  produced  from  the  OBN  receiver  functions  using  PREM  as  the  reference 
model.  Because  the  peak  on  this  R^SS  is  broad,  it  would  be  impossible  to  pick  an 
exact  value  of  R^  for  which  the  receiver  functions  stack  best.  It  is  clear,  however,  that 
the  data  stacks  best  for  values  of  Rv  less  than  the  average  1.825  value  computed  from 
PREM.  Using  the  value  of  Rv  at  the  center  of  the  peak  in  this  figure  (about  1.78)  with 
PREM  P-velocities  we  produced  the  VSSc  shown  in  Figure  17.  From  this  figure,  we 
infer  a  depth  to  the  670  km  discontinuity  of  about  675  km  with  a  fractional  velocity 
multiplier  of  about  1.01.  In  light  of  the  fact  that  the  P-velocities  inferred  from  Figure  17 
are  greater  than  those  of  the  reference  models,  the  simplest  physical  interpretation  of  a 
low  VpA/s  ratio  is  in  terms  of  higher  S  velocities  due  to  a  more  rigid  upper  mantle  than 
implied  by  PREM.  This  is  consistent  with  the  fact  that  PREM,  being  a  global  model,  is 
biased  toward  oceanic  structure:  one  may  expect  that  beneath  the  Russian  Platform 
the  mantle  is  older  and  cooler,  resulting  in  greater  rigidity.  Shearer  and  Masters, 
(1991)  give  evidence  (related  to  the  topography  of  the  670  km  discontinuity  beneath 
the  subduction  zones  along  the  northeastern  rim  of  the  Pacific  Ocean)  that  a  decrease 
in  the  temperature  of  the  upper  mantle  will  result  in  a  greater  depth  to  the  670  km 
discontinuity  and  elevated  seismic  velocities.  We  favor  the  interpretation  that  the 
depth  to  this  discontinuity  beneath  OBN  is  greater  than  670  km  and  that  the  VpA/s  ratio 
is  low  because  this  explanation  requires  a  smaller  velocity  perturbation  with  respect  to 
both  K8  and  PREM  and  it  is  consistent  with  the  Rv  value  determined  by  RvSS. 
However  it  is  clear  from  VSSc  (Figures  15  and  17)  that  for  any  reasonable  value  of  Rv 
the  average  velocity  in  the  upper  mantle  is  higher  than  that  computed  for  either  of  the 
reference  models. 

We  do  not  observe  a  200  km  discontinuity  beneath  OBN.  For  the  PREM  model,  the 
Ps  phase  from  the  200  km  discontinuity  arrives  just  after  the  P2pis  from  the  Moho 


15 


resulting  in  a  broader  peak  at  20  seconds  on  the  synthetic  VSSsph  (Figure  9)  than  on 
the  observed.  The  strong  P2pis  arrival  from  the  200  km  discontinuity  observed  at  75 
seconds  on  the  PREM  VSSsph  is  not  apparent  in  observed  VSSsph  (Figure  14). 
These  observations  lead  us  to  conclude  that  there  is  no  200  km  discontinuity  beneath 
OBN,  or  at  least  that  it  is  not  as  pronounced  as  in  PREM. 

8.  UPPER  MANTLE  DISCONTINUITIES  AT  ARU 

Our  discussion  of  the  upper  mantle  structure  observed  in  VSS  produced  for  Arti, 
Russia  (ARU)  will  draw  heavily  on  comparisons  with  the  observations  for  OBN.  The 
receiver  functions  computed  for  ARU  exhibited  a  surface  layer  response,  but  not  as 
overwhelming  as  that  observed  in  receiver  functions  produced  for  OBN.  For  the  sake 
of  comparison  of  VSS  computed  for  these  two  stations,  the  receiver  functions 
computed  for  ARU  were  low-pass  filtered  with  the  same  guassian  filter  used  for  OBN 
data. 

The  VSSsph  to  the  left  on  Figure  18  was  computed  using  the  moveout  equation  for 
the  Ps  phase;  the  one  to  the  right  was  computed  using  the  P2pis  moveout  equation. 
On  the  Ps  VSSsph.  the  observe  Ps  phases  from  the  400  km  discontinuity  (at  about  41 
s)  is  slightly  larger  in  amplitude  than  the  Ps  phase  from  the  670  km  discontinuity  (at 
about  (61  s),  which  is  similar  to  the  amplitude  relationship  observed  in  the  VSSsph 
produced  for  OBN.  This  unusually  large  amplitude  for  the  Ps  phase  from  the  400  km 
discontinuity  is  not  likely  to  be  an  artifact  of  our  method  because  this  phase  is  very 
week  to  non-exsistent  in  VSSsph  produced  for  most  other  stations  (Gurrola  et  al.. 
1992).  The  P2pis  phase  from  the  200  km  discontinuity  observed  in  the  VSSsph 
produced  from  the  PREM  synthetics  is  also  absent  in  the  P2pis  VSSsph  produced  from 
the  ARU  receiver  functions.  This  is  a  rather  surprizing  result  in  view  of  the  fact  that 
Goldstein  et  al  (1991)  observed  a  phase  that  they  interpret  as  a  reflection  from  the  200 
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km  discontinuity  in  recordings,  from  ARU,  of  nuclear  blasts  at  the  kazakh  test  site. 
Since  the  reflections  observed  by  Goldstein  et  al.  (1991)  would  be  from  hundreds  of 
km  away  (at  the  midpoint  between  the  Kazakh  test  site  and  ARU),  it  is  possible  that  the 
200  km  discontinuity  may  become  less  pronounced  closer  to  the  station.  Both  OBN 
and  ARU  are  located  on  the  Russian  Platform  whereas  Kazakh  test  site  is  over  a 
thousand  km  from  ARU  on  the  other  side  of  the  Ural  mountains.  Perhaps  this  change 
in  geological  terrain  is  the  reason  the  200  km  discontinuity  is  not  observed  in  our  VSS 
which  sample  the  mantle  within  150  km  of  the  station.  Alternatively,  the  200  km 
discontinuity  in  this  region  may  have  a  smaller  velocity  contrast  than  that  of  PREM,  and 
the  reflection  method  of  Goldstein  et  al.  (1991)  may  be  more  sensitive  to  this 
discontinuity  than  our  method. 

Figure  19  depicts  a  R\/SS  computed  for  ARU.  The  peak  in  this  figure  is  even  more 
elongated  than  on  the  RySS  computed  for  OBN,  however  the  relationship  of  a  lower 
than  PREM  (1.825)  Vp/Vs  ratio  can  be  inferred  for  this  station.  For  the  sake  of 
comparison  with  the  results  for  OBN.  we  use  the  same  VpA/s  ratio  to  compute  the 
VSSc  (Figure  20)  for  this  station  as  employed  in  computing  Figure  17.  Comparing 
Figures  17  and  20,  we  observe  that  the  peak  for  the  Ps  phase  from  the  670  km 
discontinuity  is  sharper  for  ARU  than  OBN.  More  importantly  we  observe  that  the  peak 
as  a  whole  on  Figure  20  appears  to  shifted  to  the  upper  right  corner  of  the  plot  relative 
to  the  peak  on  Figure  17  to  a  position  about  2  km  shallower  in  depth  and  no  more  than 
.01  higher  in  velocity  ratio  relative  to  PREM.  We  believe  that  this  kind  of  comparison  in 
relative  depth  and  velocity  between  two  stations  is  more  reliable  than  trying  to  pick  an 
absolute  depth  and  velocity  at  a  single  station. 
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9.  CONCLUSIONS 


Through  the  use  of  velocity  spectrum  stacks  we  can  stack  receiver  functions 
calculated  from  data  with  different  ray  parameters,  and  by  doing  so  infer  velocity 
structure  beneath  the  seismographic  station.  This  technique  can  be  used  to 
distinguish  between  a  Ps  phase  and  a  P2pis  reverberation  based  on  differences  in 
the  shapes  of  their  respective  moveout  curves.  The  method  looks  most  promising  for 
the  interpretation  of  upper  mantle  structure,  but  when  a  full  range  of  ray  parameters  is 
available,  crustal  structure  might  also  be  imaged  with  VSS.  The  shape  of  the  moveout 
curve  for  a  particular  phase  is  dependent  on  the  depth  of  the  interface  from  which  it 
originates  and  the  velocity  structure  above  the  interface.  In  order  to  compute  moveout 
for  the  curved  ray  path  of  Ps  phases  for  upper  mantle  discontinuities,  we  trace  rays 
through  a  reference  model  and  infer  a  fractional  difference  between  the  reference 
model  and  the  structure  necessary  to  satisfy  the  data.  It  should  be  clear  from  the 
examples  above  that  the  depth  of  an  interface  from  which  the  phase  of  interest 
originated  is  poorly  constrained  by  the  VSS  method  unless  assumptions  are  made 
about  the  velocity  structure  or  For  OBN  and  ARU,  we  assumed  that  the  PREM 
model  P-velocities  are  reasonably  close  to  the  truth  and  inferred  the  value  of  Rv  from  a 
velocity  ratio  spectrum  stack  (RvSS).  Applying  this  Rv  value  to  the  PREM  model 
produced  a  VSSc  from  OBN  receiver  functions  and  yielded  a  675  km  depth  for  the 
"670  km  discontinuity”  with  an  average  velocity  about  1.01  times  that  of  PREM.  The 
inferred  depth  to  the  670  km  discontinuity  beneath  ARU  is  2  km  shallower  than  that  of 
OBN  with  about  the  same  velocity  structure  as  OBN  relative  to  PREM. 

Through  the  analysis  of  VSS  produced  for  data  from  OBN  and  ARU,  we  have 
identified  upper  mantle  Ps  conversions  associated  with  the  400  km  discontinuity  that 
were  not  observable  in  the  individual  receiver  functions.  We  have  also  obtained 
evidence  that  the  200  km  discontinuity  is  not  present  or  is  weak  beneath  these 
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stations.  By  comparing  the  VSS  produced  from  PREM  synthetics  with  those  produced 
from  observed  data  we  conclude  that  the  greater  velocity  contrast  across  the  400  km 
discontinuity  proposed  in  the  K8  model  (Given  and  Helmberger,  1980)  for  upper 
mantle  structure  beneath  northwestern  Eurasia  is  more  appropriate  for  this  station  than 
that  of  PREM. 
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Figure  1 .  Ray  paths  for  the  Ps  phase  relative  to  the  P  phase  for  a  layer  over  a  half 
space.  Tp  and  Ts  are  the  travel  times  of  the  P  and  S  phases  with  the  same  ray 
parameter  through  the  layer  respectively.  Th  is  the  travel  time  differential  in  the  half 
space  for  the  two  rays  assuming  a  planar  wave  front. 
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Figure  2.  Seismic  section  of  synthetic  receiver  functions  computed  for  a  range  of  ray 
parameter  from  0.040  to  0.080  by  raytracing  through  a  model  with  a  40  km  thick  layer 
(Vp=6.0  km/s,  Vs=3.5  km/s)  over  a  half  space  (Vp=8.0  km/s,  Vs=4.6  km/s). 
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Figure  3.  Ray  paths  for  the  P2pis  phase  relative  to  the  P  phase  for  a  layer  over  a  half 
space.  Tp  and  Ts  are  the  travel  times  of  the  P  and  S  phases  with  same  ray  parameter 
through  the  layer  respectively.  Th  is  the  travel  time  differential  in  the  half  space  for  the 
two  rays  assuming  a  planar  wave  front. 
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Figure  4.  Velocity  spectrum  stacks  produced  from  the  synthetic  receiver  functions 
depicted  in  Figure  2.  Ps  stacks  are  shown  on  the  left  and  P2pis  stacks  are  on  the 
right. 
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Figure  5.  Stacks  of  the  synthetic  receiver  functions  depicted  in  Figure  2.  On  top  is  a 
straight  stack  with  no  time  correction  applied.  The  bottom  three  receiver  function 
stacks  are  computed  after  applying  the  appropriate  normal  moveout  correction  for  Ps. 
P2pis,  and  Pip2s  respectively  (from  top  to  bottom). 
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Ps  Moveout  Curves 


Figure  6.  Moveout  curves  at  various  depths  from  25  km  to  700  km  for  the  Ps  phase 
computed  for  the  PREM  velocity  structure.  The  solid  lines  are  moveout  curves 
computed  for  a  single  layer  with  equivalent  average  slowness  as  PREM  to  the 
respective  depths.  The  dashed  line  is  the  moveout  computed  for  a  curved  ray  path  by 
ray  tracing  through  PREM.  The  depths  associated  with  each  curve  are  printed  on  the 
left. 
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P2pis  Moveout  Curves 


Figure  7.  Moveout  curves  at  various  depths  from  25  km  to  700  km  for  the  P2pis  phase 
computed  for  the  PREM  velocity  structure.  The  solid  lines  are  moveout  curves 
computed  for  a  single  layer  with  equivalent  average  slowness  as  PREM  to  the 
respective  depths.  The  dashed  line  is  the  moveout  computed  for  a  curved  ray  path  by 
ray  tracing  through  PREM.  The  depths  associated  with  each  curve  are  printed  on  the 
left. 
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Figure  8.  A  representative  display  of  the  80  low  pass  filtered  (corner  frequency  of  0.3 
Hz)  PREM  synthetic  receiver  functions  used  to  compute  the  VSS  in  later  figures.  The 
moveout  curve  associated  with  Moho  and  upper  mantle  discontinuities  are  labeled. 
The  dashed  lines  are  the  true  moveout  curves  produced  by  raytracing  through  PREM. 
The  solid  lines  are  the  moveout  curves  computed  by  replacing  PREM  with  a  layer  over 
a  that  would  result  in  a  vertical  travel  time  equivalent  to  that  of  PREM. 
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Figure  10.  The  VSSc  produced  from  the  PREM  synthetics  depicted  in  Figure  8.  The 
VSSc  on  the  top  was  computed  using  PREM  P  and  S  velocities  as  the  reference 
model.  The  reference  model  used  to  produce  the  two  VSSc  on  the  bottom  used  P 
velocities  from  PREM  (left)  and  K8  (right)  and  a  value  of  1 .825  for  Rv. 
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Figure  11.  Velocity  ratio  spectrum  stacks  (RySS)  computed  using  PREM  (left)  and  K8 
(right)  reference  models  for  the  synthetic  receiver  functions  generated  for  the  PREM 
model. 
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Figure  12.  A  representative  sampling  of  the  113  low  pass  filtered  individual  receiver 
functions  computed  for  OBN.  The  only  clear  arrivals  without  stacking  are  the  Ps 
conversion  from  the  Moho  •seconds  and,  in  a  few  receiver  functions  the  P2pis  and 
Pip2s  from  the  Moho  at  2.^  and  25  seconds  respectively. 
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Figure  13.  Synthetic  (solid  line)  and  observed  (dashed  line)  stacked  receiver 
functions  computed  for  the  IRIS/IDA  Seismographic  station  at  Obninsk.  The  synthetics 
pertain  to  the  crustal  structure  models  depicted  on  the  left  of  the  respective  receiver 
functions.  The  low  pass  filtered  response  (stack  of  the  receiver  functions  given  in 
Figure  12)  is  given  in  the  top  frame;  the  broad  band  results  are  in  the  bottom  frame. 
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Figure  14.  Velocity  spectrum  stacks  produced  from  the  OBN  receiver  functions  shown 
on  Figure  8.  The  plot  on  the  left  was  computed  for  Ps  phases  while  the  one  on  the 
right  was  computed  using  the  P2pis  moveout  curve. 


Figure  16.  Velocity  ratio  spectrum  stack  (/?, 
computed  using  PREM  as  the  reference  model 
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Figure  18.  Velocity  spectrum  stacks  produced  from  receiver  functions  computed  from 
seismograms  recorded  at  ARU.  The  plot  on  the  left  was  computed  for  Ps  phases 
while  the  one  on  the  right  was  computed  using  the  P2pis  moveout  curve. 
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Figure  19.  Velocity  ratio 
computed  using  PREM  as 


depth  (km) 


Figure  20.  VSSc  computed  from  ARU  receiver  functions  similar  to  that  shown  in 
Figure  19,  except  in  this  case,  the  V£ 
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